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THEORY  O’7  THIN  WING  IN  A SUPERSONIC  ^LOW  WITH  CONSIDERATION  OP  THE 
NON-EQUILIBRIUM  STATE  OP  EXCITATION  0^  OSCILLATING  DEGREES  OP 
FREEDOM 

Ye,  P,  Aksenov,  Yu,  M.  Gripor'yev 

Tn  rapid  pas  Plows  we  can  observe  chenical  anti  thermodynamic 
non-enuilibrium  state.  It  is  known  that  after  the  action  of  per- 
turbation, the  propressive  and  rotarv  components  of  enerpv  rapidly 
assume  their  eouilibrium  values,  and  the  oscillatinr  component  of 
enerry  achieves  its  eouilibrium  value  many  times  slower;  and  the 
relaxation  time,  as  we  call  this  time  interval,  can  prove  to  be 
substantial.  This  permits  us  to  take  the  followinp  scheme  of  exam- 
ininp  the  non-enuilibrium  state  in  our  problem. 

Tn  a shock  wave  there  occurs  an  initial  excitation  of  the 
propressive  and  rotarv  deprees  of  freedom  under  conditions  of  a 
^rosen  state  of  the  oscillatinp  depree3  of  freedom. 

After  the  shock  wave  we  have  non-enuilibrium  excitation  of 
oscillations  under  those  conditions  where  there  is  a place  for 
the  eouilibrium  state  between  the  propressive  and  rotary  deprees 
or  freedom. 
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It  is  proposed  that  total  energy  Ea  of  active  (i.e.  progres- 
sive and  rotary)  degrees  of  freedom,  anvwhere  in  the  gas,  has  an 
enuilibrium  value  of  Ei=C‘,,7’,  and  oscillating  energy  Ev  satisfies 
relationship  [2]: 

Hr  (?,-£,)  , (D 

v.'here  Cv  - enuilibrium  value  of  oscillating  energy,  and  - 

relaxation  time. 

The  purpose  of  the  present  work  is  the  calculation  of  the  non- 
enuilibrium  state  in  excitation  of  oscillating  degrees  of  freedom 


in  the  problem  of  the  steady 


f lov.'-around  of  a thin  wing  of 


infinite  span  by  a supersonic  gas  flow. 

A system  was  taken  as  the  initial  svstem  of  eouations,  con- 
taining the  usual  enuations  of  r-as  dynamics  with  the  addition  of 
a relaxation  eouation  (1). 

As  applied  to  our  problem,  this  system  has  the  form: 

v»-3F  + tV“3 J-~TTx  ' 

dvy  , * dvy  l dp 

T~d~y  • 

'&+*£+»(fc-+4iL)-<>-  ' (2) 

«■  Iff  (d"  r+ £. ) + (<$’  r+E.)  - -f  (£  +$) . 

= 0,»  — 

P-RpT. 

In  these  equations,  t>„  vy  are  the  components  of  macroscopic 
gas  rate,  p - density,  p - pressure,  T - absolute  temperature, 

C*'*  - specific  heat  with  a constant  volume,  relating  to  the  active 
degrees  of  freedom. 
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The  desired  hydrodynamic  elenents  should  satisfy  the  boundary 
conditions  on  an  unknown  line  of  strong  disturbance  and  the  flow- 
around  conditions. 

Let  U3  propose  an  approximate  analytic  solution  to  this  prob- 
lem which  rives  prepared  formulas  for  computing  the  desired  func- 
tions in  any  point  of  the  profile.  With  the  solution  to  the  prob- 
lem in  the  first  approximation,  the  acnuired  formulas  are  similar 
to  the  known  Akkeret  formulas. 


1.  Linearized  equations  and  boundary  conditions  of  the  problem. 

Let  us  examine  a thin  slightlv-curved  profile  with  sharp 
edges  at  small  attack  angles. 

We  introduce  into  the  system  a coordinate,  the  beginning  of 
which  we  place  in  the  front  part  of  the  wing;  axis  x we  set  along 
the  inflowing  flow,  and  axis  y above  and  perpendicular  to  axis  x. 

Since  the  thin  profile  barely  disturbs  the  inflowing  flow, 
the  desired  hydrodynamic  elements  can  be  given  in  the  form: 

(1.1) 

P — Pt  +P 

r=  r,  + r, 

Et  = ( r,)  + £/, 

p = Pi  + p\ 


where  p^,  V^,  p^,  and  - values  of  hydrodynamic  elements  in  the 
undisturbed  flow  - constant  values. 

v*.  •/.  ? > T* , £»',  p'  and  those  produced  along  the  coordinates  - 

smalls  of  the  first  order. 
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Linearizing  system  (2),  i.e.  disregarding  the  terms  above 
of  the  first  order  of  smallness,  we  obtain  system: 


. „ \ */>' 

1 dx  p,  1£T  • 


(1.2) 


v */  = 1 y 

K'  dx  p,  * 


v',  — jg-  = u>o  (Cv  ’ V — £/), 


/p'-/?(p1r+rl'p')i 


where 


*-(#),  • 


For  findinp  the  functions  of  vt',Vy,p\  T,ES'.  p'  which  interest 
us  by  integrating  system  (1.2),  we  must  write  down  the  boundary 
conditions  for  these  functions  which  interest  us. 

Let  U3  bepin  from  conditions  on  the  surface  of  the  nonremov- 
able discontinuity.  We  a free  to  provide  by  subscript  1 the  values 
of  hydrodynamic  elements  up  to  the  nonremovable  discontinuity/,  and 
by  index  2 the  values  of  hydrodvnanic  values  after  teh  surface  of 
nonremovable  discontinuity.  Then  the  momentum  theorem,  the  mass 
conservation  lav/,  and  the  law  of  conservation  of  enerpy  for  pas 
masses  passinp  through  the  shock  wave,  are  written  down  in  the 
form  of  a relationship: 

— Pi  rH„(V„a  — V,.i)  = (Pi  ~ Pi)COS  («,*), 


— Pi  V,„i  (V,,i — Vy.,)  = (Pa~  Pi)  COS  (ft,  y), 

Pj  I \\,t  cos  (ft,  x)  4-  Vy„  cos  (rt,y)|  - p,  vn„, 

-f,v..,(iL±ik+E,- *■;,  -g,)- 

= Pi  I Wj  cos  («,*)[  + ryilcos  (ft,  y)  1 — p,  v‘n„,' 


(1.3) 


\ 


-4 


mm 


where  V ^ ■VK  |Co*(n,  x)|+ Vy>i  Cos  (n,  y)  projection  of  speed  \T^  of  an 
undisturbed  flow  on  a normal  n of  the  surface  of  nonremovable  dis- 
continuity . 

Since  the  x-axis  is  directed  according  to  an  undisturbed 

flow,  vXji  * Vi»  vyti  a °*  If*  we  designate  through  <p  the  angle  be- 
tween the  tangent  of  breakdown  and  the  x-axis,  then  Cos  (n,  x)-Sin<p, 
Cos  (n,  y)-—  Cosf<p  , and  consequently  v'n.i  = k,  sin  ? . 

Then  from  (1.3)  we  find: 

l 


V*~  * vr**  = 7" fy(Pi—Pt). 


H W»in«t -(*-/,,) 


Pt’  — P,’ 


Let  us  propose  ? = « + rAe sin « = — = ~ 

g,  M | 


(1.4) 

(1.5) 


a^  - speed  of  sound  in  a disturbed  flow. 


Let  us  consider  that  A«p  - value  of  the  first  order  of  small- 
ness. Then  with  accuracy  up  to  smalls  of  the  first  order,  we  have 


sin*  9 = sin*  («  + A?)  = — !—  A*. 


ctg  f = VXi,*—  i — M*  A<p. 


(1.6) 


Taking  into  account  (illegible)  and  (1.6)  from  the  relation- 
ships (1.4),  we  obtain  for  the  desired  functions  the  following 
conditions  on  the  line  of  noremovable  discontinuity: 


(1.7) 

_ Vas- 1 , 
h V p ' 

(1.8) 

(1.9) 

* = grr!>. 

(1.10) 

> j 
1 

\ 
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Moreover,  considering  the  frozen  state  of  oscillating  derrees  of 
freedom  with  the  pnssinr  of  gas  through  the  shock  wave,  we  have 
one  condition  on  the  line  of  discontinuity 

Ey  =0.  (1»11) 


Let  u 3 direct  our  attention  now  to  the  condition  on  the  pro- 
file (to  the  flow-around  condition). 

We  have 


(1.12) 

where  p is  the  angle  of  inclination  of  the  tanrent  to  the  profile 
to  the  x-axls. 

Let  us  linearize  this  condition  on  the  profile.  Let  the  eoua- 
tion  of  the  profile  be 


>-«<*>•  (1.13) 

Since  the  profile  is  thin  and  mildly  curved,  it  means  that  value 
tRP-C'(x)  can  be  considered  snail  of  the  first  order.  Assuming  in 
(1.12)  vy  = vy\  sb  V,  + t>/  and  disrerardinr  the  value  of  the  small- 
ness of  the  second  order,  we  obtain  condition  on  the  profile  in 
the  form 


v- *«'(*>.  (i.no 

Thus,  we  nu3t  integrate  the  system  of  equations  (1.2)  with  boundary 
conditions  (1.7 )-(l.  11 ) and  (l.HO. 

2,  deneral  solution  of  the  linearized  system  of  equations. 

Let  U3  reduce  the  system  of  equations  (1.2)  to  one  enuation 
relative  to  function  p*.  With  this  purpose,  we  will  differentiate 
the  fifth  enuation  of  system  (1.2)  twice  along  the  x-axis.  We  ob- 
tain8 • (2.1) 
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^r-orr  teh  fourth  eouation  of  system  (1,2)  we  find 

0£t'  Pl  ( dt-,'  , dv \ „„>  dT’ 

Sx~~  7 TPT'  dx  v “31“  ' 

Differentiating  (2.2)  once  with  respect  to  x,  we  obtain 


Sir 


a_  *jv \ _ C(.)£ZL 

Pt  Pt  ' d.xlfy ) v dr3 


(2.2) 


(2.3) 


The  third  equation  of  system  (1.2)  with  the  use  of  the  latter  pro- 
vides 

dT"  ___  I dp’  T,  / jf,'  dt'y'  . 

dr  ~ /?  p,  dx  ■*"  T7  ' Sx  ’ dy  ’ ‘ (2  k) 


(2.4) 


Having  differentiated  (2.4)  once  with  respect  to  x,  we  find 

JLL- __*_**£'  i 7-.  , d>rv’\  (r> 

Sir-  Rfi  dxz  -i  Vi  (-55T  + -3 

Let  us  differentiate  the  first  equation  of  system  (1.2)  with  re- 
spect to  x,  and  the  second  with  respect  to  y.  We  find 


d * f,' 

1 

d'-p’ 

(2.6) 

3J»  ~ 

Pi  V'i 

d-«»  ’ 

d>  Vy’  _ 

1 

d*  p’ 

(0  7\ 

dxdy  ~ 

Pi  V'i 

dy5  ’ 

?lacinr  (2.5),  (2.6),  and  (2.7)  in  (2.3),  we  obtain 


) 

in 

(2. 

T\ 

4-  - 

_£hJ 

Pi 

i 

Pi"  / 

_£i 

- 4- 

Pi 

i 

Pi*  / 

(2.6) 


Differentiating  (2.3)  once  with  respect  to  x and  placing  the  ob- 
tained equation,  along  with  (2.5)  and  (2.3)  in  (2.1),  using  (2.6) 
and  (2.7),  we  obtain  for  function  p’,  after  the  transformations  of 
coefficients,  and  equation  with  partial  derivatives  of  the  third 
order  in  the  form: 

*"a**^r  -ko£§?  + b's£--sfr  = 0'  (2.9) 

where  the  following  designations  are  introduced: 


(see  following  page) 
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*.-*_*iS 


<•> 


r + <•»+  d*'  ' 


-(•» 


(J*  = 


R T,  (Af  + c‘*‘) 


l'.»- l = »,•-!  , 


c<;,+c,;> 


*r,  <# + c["+  c<;,‘) 


■—T  y,»—  1 * M,1  — 1 . 


(*) 


Equation  (2.9)  replaces,  In  the  non-equilibrium  cane,  the  normally 
used  wave  equation.  Tills  tie  with  the  classic  equation  is  quite 
apparent.  When  the  relaxation  time  approaches  null,  from  (2.9), 
as  in  the  particular  ca3e,  an  equation  follows  °or  an  equilibrium 
period.  In  the  case  whereiko  -*  oo(inf initelv  larper  time  o*'  relaxa*- 
tion),  we  obtain  a special  c*3e  of  the  classic  equation  with  the 
"nch  number  of  the  "frozen"  period,  "or  some  non-null  finite  value 
of  relaxation  time  the  characteristics  of  equation  (2.9)  with  the 
form 


play  the  sane  role  as  the  characteristics  which  are  well  known  for 
equilibrium  periods. 

We  will  search  for  a solution  to  the  obtained  equation  (2.9) 

In  the  form  of  an  exponential  law,  which  is  natural  for  the  relax- 
ation proce S3. 

It  Is  easy  to  see  that  the  functions  with  the  form  e«*  — |»y and  e«** 
— >*v,  where  a mA«*.  fl  arJL  P“  are  complex  con.lurate  numbers,  are  solu- 
tions to  equation  (2.9)  with 


(2.10) 


.’hen  function 


■-&  u-  P*y 


(2.11) 


where  \ and  P .ire  the  random  constants,  and  will  also  be  the  solu- 
tion to  equation  (2.9). 


e rropose 


« = p+q  l,  P = J + X / 


'hen  (2,11)  can  be  written  in  the  form 


p'  = t px  Jy  |Ccos(?jc  — Xy)  -f  Dsin  ( qx  — X y)J, 


(2.12) 


(2.13) 


where  0=\+h,  D»1(A-n).  Substituting  (2.12)  in  (2.10)  and  separating 


the  real  and  irapinary  rarts,  we  obtain 

**-X«  = F, 
2 »X  = G, 


v'here 


f-  ((»» + »,«V)  (;»-?»)  - 2 fc, (■ 

(l  + + *0*  h 

1 q (P*  — q')  4-  (»»  + kna'p)2pq\  k„q 

(1  +*pP)>  + *„>?’ 

n _ 1^«» «(j*-g»)4(»»  + *o«»;>)  2 />g)(i  +*,/>) 

(1  + *»  pY  + kjf 
(jP -«»)- 2 *.«W1  **. 


J"T,  (2.14)  we  find 


(i  + M)*+*-y 


*,  = -5-(/r+K7rT5f).  x*  = -j-(-/r+  K/^  + g*). 


(2.14) 


(2.15) 


(2.16) 


Havinr  por  "unction  p'  the  expression  (2.13),  we  find,  proceedinr 
bv  opposite  means,  the  remaininr  hydrodynamic  elements : v/.  Vj'.  V,  p',  Ev'. 
Intei-ratinr  the  "irst  equation  of  system  (1.2),  we  find 


*«'  = — 77vV  e?X  ly|CcOS(^X-Xy)  + 
-f-  D sin  (qx  — Xy  ) ] -(-  C,  (y). 


(2.17) 


where  C^(v)  is  the  arbitrary  function, 
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Prom  the  second  equation  of  system  (1.2),  placing,  instead  of 
l>' , its  equation  (2.13),  we  obtain 

dJto  =~7rW  ePX~ty  I(-*C-XD)Cos(?JC-Xy)  + ^ ^ 

+ (XC  — ID)  sin  (qx — X y )J  . 

Integrating  (2,18)  with  respect  to  x,  we  find 

vr  = WV,(p‘  + q ) Q*—*9)—D  (Xy+ »;>)]  sin(vJf-).y)+  (2.19) 

+ \C(—hp  — \q)  + D(ty-f  Ip)]  cos  (qx  — Xy)}  + C,(y), 

where  02(y)  - arbitrary  function. 

^rom  the  third  equation  of  system  (1.2),  placing,  instread  of  «/ 

y/ , their  expressions  (2.17)  and  (2.19),  we  integrate  with  re- 
spect to  x and  find 


P1  - -^77— I Ccos  (qx-ly)  + D sin  (<7X-Xy))  -$;C,'(y)jc  + 

+ + Sin  (qX  ~ W + 

+ ( Cf-Dg ) cos  (qx  - Xy)J  + C,  (y), 


g = 2pq  F J-  (q'  — p*)  G, 
f=  (p* -q*)F+2pqG, 


and  C^(v)  - arbitrary  function. 

From  the  last  equation  in  system  (1.2)  we  find 


r = — ft),  p,»  * <fPX  iy  I C COS  (qx  - Xy)  + 
+ Dsin  (qx  — Xy)]  -f  C2'  (y)-jf- 

■ \(Gg  + Df)  sin  (qx  — Xy) 


r ?px  ly 


Pi  +,?»)* 


-I-  (Cf  - Dg)  COS  (qx  - Xy)j  - J±.  C,  (y). 


(2.21) 


(2.22) 


'ron  fourth  equation  of  system  (1.2),  usinp  the  fifth,  we  find 

£-'  - [^^r-  - ,7^;]  ^ " *r ' sl"  <♦*- 


.<«>  T 
«v  T> 


+ : 


MO  + ( Cp  + Dq)  t Q3  (qx  - Xy))  - ■ 

Px  \ <px“*y 

^Tp,  j -p » + qt  l(C«  + Du)  sin  (qx  — Xy)  + 
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ja-*y 

+ (Cm  - Dm)  cos  ( fx  - Xy)|  - |(Cf  + 0/)  X 

, f r« 

X*la(f*-Xy)+  (C/-Df)co»(f.*-Xy)l+  ^-»^— + 


c°,)  7i  p, 

+ -VTx  + lfc 


■)C* 


'(y)  + 


4-  .*r - *IL*px  ->y  (Cco«  (vx  - Xy)  + £>iia  (qx  - Xy)j  - 


where 


J«>r 

C*  (y)  ( 


(2.23) 


pO, 

n = pF+qG.  (2.2b) 

p't  v,',  v,*,  p'  T'  £L' 

^unctions  Riven  in  formulas  (2,13),  (2.17),  (2.19), 

(2.20),  and  (2.23),  satisfy  the  five  equations  of  system  (1.2). 
Satisfying  the  fifth  equation  of  system  (1.2),  we  obtain  one  rela- 
tionship betv/een  the  random  elements  entering  into  the  equation. 

So,  in  the  general  solution  to  system  (1.2)  arbitrary  functions 
G^(y),  C2(y),  and  C^(y)  enter,  as  do  the  random  constants  C,  D,  n, 
and  n,  which  we  must  find  from  the  boundary  conditions. 

3.  Solution  to  the  problem  in  linear  approximation  (search  for  ar- 
bitrary functions). 

We  find  arbitrary  functions  which  entered  into  the  general  so- 
lution. "or  this,  we  have  six  boundary  conditions  (1.7 )-(l.ll) , and 
(1.14). 

Conditions  (1.7 )— (1. 11)  on  the  line  of  nonremovable  discontinuity 
will  be  written  on  characteristic  x - ay  * 0,  and  condition  (1.14) 
will  be  on  the  y-axis  y ■ 0, 

Satisfying  condition  (1,7)*  we  see  that 

C,(i)-0.  (3.1) 
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nnd  Prom  condition  (1,8)  wo  find 

■jlP*-*) 

o. (t) - *hr+w 11 |C <■** + ^ ~ * ~ 1 + D (‘* " 

- lp\\  CO*  -j  («f  - *)  + 1C  (l#  - If)  + D(aq*  + ap*  — lp  — 
-Xf)l«»n 

Satlsfyinr  condition  (1.9),  we  rind 

— C,(y)-  (Tj-7r)#y<P*~')  |Cco*y(«f— 8)  -f-  £>»lny  (aq- 
y <?*-»> 

— Mil  V’.V/'./rv7)'  1 £ t -Df)  sln(flv~My  \-(Cf-Dg)  Cos  aq - 


" Myl  — IVfflW1)1 1 CM  4 DN)  cot  (aq  — K)y—  ( CN 
— DM)  sin  (aq  — 8)  y | a y , 


where 


l 

M — a{p»  ff  </*)  (ap  - 28)  4 pF  -4-  qG, 
N « «i(/>*4  v*)(av  — 2M  —<tF  + pO. 


Condition  (1.11)  loTvIn  to  count  Ion 


where 


e [^C  + BD)  cot  ( qx  - X i- j + [AD  - BC)  *in  (qx  - 


!-/>(.•-  I)( 


(.  f*L 

P\ 

\ «■»*  K 

+ 7rZv: 

O' 


M 

| C»  n ..  f 

1 7; 


(•) 

V^Vr-2aX, 


which  will  be  rmtinfled  1C  we  net 


— 0.  B — 0. 

Uni  nr  condition  (1.1M),  we  obtain 

-j-Pi  ^i*?'  (-*“)  *(C co*  | Dtlnqx). 


( 3.0 

(3.3) 

( 3 JO 

( 3 . '<) 

( 3 . (> ) 

(3.7) 

( 3.8) 


1C 


(3.9) 


deplacinr  in  equation  (3.8)  qx  by  qx-Ay,  we  iy.nd 

-p*  + piy 

c cos  (qx  — Ay)  4-  £>tin  (qx  — Xy)  — -Cl-L. e 


<'(*- 


then 


(3.30) 


The  obtained  expression  for  p’  is  a renerallzat ion  on  the  cane 
of  movement  of  ran  under  conditions  of  the  absence  of  thermodynamic 
equilibrium  of  the  known  ARkeret  formula  for  the  thin  winr  and  con- 
curs with  *.-0  W *0-“  . 

Alonr  the  3a me  lines,  in  the  case  *»-*•()  from  (2.1)  it  follows 
that  p*— ■*<**  and,  conseouentlv, 


k_ 

* 


(3.11) 


py  (---I 

f>o  that  multiple  e * * , considerinr  the  non-equilibri- 

um state  in  this  case  turns  to  a unit,  we  obtain  the  Akkeret  for- 
mula, In  the  case  *,-*■« 

from  (3.10)  we  obtain  the  Akkeret  formula  with  a ".ach  number  of 
the  "frozen”  period. 

Knowing  n*,  we  can  compute  the  coefficients  of  lift  force  and 
dra.r  resistance  and  determine  in  this  manner  the  same  contribution 
to  the  values  of  these  coefficients,  which  introduces  the  non-enui- 
libr.ium  state  of  excitation  of  oscillatlnr  derrees  of  freedom. 
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